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INTRODUCTION 


Two significant facets of Apollo oxygen tank 
operation are the occurrence of thermal stratification at 
high fluid densities and the tendency toward high heater 
temperatures at low fluid densities. Some insight into the 
nature of these phenomena can be gained, respectively, by 
consideration of a two fluid stratification model and a 
conductive heat transfer model. The alternative of convec- 
tive blowdown tank operation at low densities is briefly 
examined. 


STRATIFICATION 


Thermal stratification is significant in the 
Apollo oxygen tanks because of the greatly reduced mixing 
in the low gravity field and the low thermal conductivity 
of the oxygen. In the stratified tank adjacent phases of 
different temperature coexist stably. The less dense phase 
exerts a pressurizing influence on the relatively incompres- 
sible dense phase. The effects of stratification are for 
the most part limited to high density conditions; for 
example, pressure decays have not been observed on Apollo 

flights for densities of less than M2 lb/ft 3 (60% tank 
quantity) . The phenomena of pressure decay relates to the 
fact that at the higher densities the energy input neces- 
sary to achieve a given pressure rise is measurably less 
with stratification than without. Hence the tank can main- 
tain operating pressure despite a small energy deficit; 
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induced mixing, as caused by a vehicle maneuver, then 
results in a sudden pressure drop or decay. 

By conservation of energy the relation between 
energy input, mass outflow, and pressure rise rate for a 
homogeneous fluid is: 
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and Q, to, V, p, u, and h are the heat input, mass outflow, 
volume, pressure, internal energy, and enthalpy, respec- 
tively. The coefficient C accounts for the change in tank 
volume with increasing pressure. Equations (1) - (4) 
yield the pressure rise rate when heat is added uniformly, 
that is, stratification does not occur. 

Over the entire pressure cycle the energy input 
is related solely to the fluid's enthalpy-density relation- 
ship 0 and the rate of mass outflow from the tank to. That 
is, there is no net pressure change so Equation (1) becomes 


Q = (00 


(5) 


Figure 1 shows the energy input data for Apollo 14. The 
solid curves were computed by Equation (5) for a mean tank 
pressure of 900 psi using the average flow rates and heat 
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leaks indicated. Tank quantity rather than density is 
used as the abscissa, with 100% quantity equal to 

3 

69.5 lb/ft . The spread in the data reflects the variation 
of the flow rates about the mean value. 

The time to pressurize the tank within its 
operating band offers a direct measure of the stratifica- 
tion occurring. Figure 2 portrays this data for Apollos 12 
and 14. The spread of the data points below 60% for the 
two missions is traceable to two effects.* Firstly, 
because of the tank arrangement on Apollo 14, Tank 3 mass 
outflow during its heater operation was substantially 
increased by leakage of the Tank 2 check valve. In effect 
Tank 3 pressurized Tank 2 during each cycle. Secondly, for 
thermal reasons only two of the three Tank 3 heater elements 
were employed below 41% (with two exceptions). With this 
in mind the pressurization times for Apollos 12 and 14 
are essentially identical even though Apollo 12 utilized 
fans for periodic mixing (twice daily). That is, the fan 
operations did not appear to noticeably reduce stratifica- 
tion effects on pressurization. 

In Figure 3 the Apollo 14 data are plotted with 
the pressurization time predicted by Equations (1) - (4) 

indicated by the solid lines. The flowrates employed are 
given in Figure 4. The flow values are taken from the 
quantity guage readings, corrected in the case of Tank 3 
for the flow necessary to pressurize Tank 2. As expected, 
the agreement between data and prediction is relatively 
good below 60%, but not above. 


*Note: On Apollo 14 all data above 60% is from 

Tanks 1 and 2 and below 60% is from Tank 3; whereas on 
Apollo 12 Tanks 1 and 2 were operated concurrently through- 
out (there was no third tank) . 
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The basic elements of the effect of stratifica- 
tion on pressurization are the localized storage of energy 
in the vicinity of the heater and the nature of the 
enthalpy-density relationship as reflected in 0. Following 
heater activation the fluid adjacent to the heater 
increases in temperature, and hence expands. This expan- 
sion mechanically pressurizes the remainder of the fluid. 
The pressurization of the bulk fluid by the hot "bubble" 
is given by 
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where the tank heat leak and outflow have been neglected. 
The bulk's volume change is equal and opposite to that of 
the bubble, whose mass is assumed constant: 
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The stratified pressure response of the system, Ap, to a 
localized heat input, AQ, is therefore given by: 
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with subscripts s and B referring to the stratified bubble 

and bulk fluid, respectively. The 0 curve passes through a 

. . 3 

minimum at a density of ^26 lb/ft . If the bulk density 

is greater than this, the 0 and p variations accompanying 

the creation of a less dense phase (film) adjacent to the 

heater combine to cause a larger pressure rise for a given 
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3 

heat input. Below ^42 lb/ft density (60%) the 0 decrease 

3 

with expansion falls off and below 26 lb/ft the 0 varia- 
tion counteracts the effect of density decrease on pressure 
rise. 


A rough indication of the extent of the stratifi- 
cation process can be gained by a two fluid calculation. 

The model is schematically depicted in Figure 5. Initially 
all heater energy Q HTR is assumed confined to the strati- 
fied volume adjacent to the heater. The energy is accom- 
modated by fluid outflow m^ from this volume. The bulk 

fluid receives this stratified outflow as well as heat leak 
Q H lk f r °m the environment and supplies tank outflow m Q to 

the spacecraft systems. Equations (1) - (4) describe the 
pressure rise in each phase. When the stratified fluid 
density has diminished to that point where its 0 value has 
passed through the minimum and risen to the bulk value, 
the heater energy is assumed to be transferred uniformly 
to the bulk phase. Of course, other phenomena take place, 
such as stratified volume growth, mass transfer from the 
bulk, and energy sharing with the bulk. Nevertheless, 
the simplified model exhibits the basic attributes of the 
observed pressure behavior as shown by the dashed line in 
Figure 3. The calculation results are for a stratified 

3 

volume of 0.0119 ft (1/4% of tank volume) or a radial 
distance of 1/8 inch. This is on the order of boundary 
layer dimensions, and explains the insensitivity of 
pressurization times to overall tank flow condition as 
evidenced by the close grouping of Apollo 12 and 14 data 
for a variety of spacecraft inertial conditions. 

Some interesting aspects of stratification are 
shown in Figure 6. Here fan activation causes the expected 
pressure decay. The divergence from uniform heating condi- 
tions of the pressure profile prior to fan operation 
reflects stratification. The pressure behavior during the 
heater cycle following fan operation is significant, 
though. As expected, in the well mixed tank the pressure 
behavior follows a uniform heating path. However, the 
recovery from the low pressure point occurs at a rate 
identical to the pressure rise rate prior to fan operation. 
Hence, even in this clearly active fluid motion environment 
some stratification of the same type as that observed in 
the quiescent period prior to fan operation occurs. 
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CONDUCTIVE HEAT TRANSFER MODEL 


Since the gravitational forces exerted on the 
vehicle during steady flight are quite small, ranging from 

_ 0 _g 

roughly 7 x 10 g in attitude hold to 3 x 10 g during 
3 rph passive thermal control, the buoyant motion of the 
fluid during heat transfer periods is correspondingly 
reduced. In the limit, heat is transferred into the fluid 
solely by conduction. This constitutes a worst case in 
terms of maintaining a low heater temperature while trans- 
ferring the heat required to sustain tank operating pres- 
sure. Therefore, a conduction model was developed to 
predict worst case (zero gravity) heater temperatures 
(Reference 1). 

The simplified model is schematically depicted in 
Figure 7. It assumes spherical symmetry (as the tank) and 
adapts the cylindrical heater parameters to this configura- 
tion. The spherical heater area is set equal to its 
cylindrical counterpart, while its internal volume is 
adjusted to match the true cylindrical volume. Twelve 
fluid nodes of equal thickness are dispersed between the 
heater and tank wall; the fluid within the heater is 
represented by one node; and the heater, tank wall and SM 
environment are represented by one node each. 

The CINDA-3G numerical heat transfer program was 
employed. The logic of this program is shown in Figure 8 
via its electrical analogue where each node with its own 
capacitance is linked to adjacent nodes by resistive cir- 
cuits. The radiative transfer between the heater and the 
tank wall is included as is the heat leak to the tank from 
the surrounding environment. Tank pressure is held constant 
throughout each calculation, and the nodal mass, specific 
heat, and thermal conductivity are adjusted at each time 
step in accord with the nodal temperatures resulting from 
the conductive heat transfer. The oxygen properties as 
correlated by R. B. Stewart (Reference 2) are employed for 
density and specific heat, while the thermal conductivity 
is taken from a North American Rockwell report (Reference 3). 

The model is admittedly crude. It does, however, 
furnish some interesting results. Firstly, the predicted 
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temperatures agree well with those for steady attitude 
hold heater cycles. Figure 9 illustrates this. Here the 
temperature profile and maximum temperature indicated by 
the heater temperature sensor are closely matched by the 
model. (One qualification to be made is that the model 
assumes a uniform heater temperature while the data is 
from a thermocouple linked to one point on the heater. 

The excellent correlation between the two is taken to mean 
that the temperature sensor reflects an effective mean 
heater temperature from a heat transfer standpoint.) It 
appears that some mixing occurred at the beginning of the 
heating cycle but after this delay the temperature rise and 
fall off after heater deactivation are as predicted. Thus, 
the g forces exerted during attitude hold seem to produce 
little convective mixing, although the more rapid return 
of the heater to its original temperature indicates some 
convection. 

It is interesting to look next at the ability of 
spacecraft rotation to enhance heat transfer.. Figures 10 
and 11 show the heater temperature profiles for rotation 
rates of 1 rph and 3 rph respectively. The ratio of the 
maximum temperature rise achieved in each case to the 
maximum temperature rise predicted by the conductive model 
offers a rough measure of the convective heat transfer 
coefficient. The increased rate of rotation clearly 
produces improved heat transfer. 

Another form of convective motion is shown in 
Figure 12. Here the spacecraft has been brought to rest 
from a 3 rph passive thermal control rotation. The low 
temperatures indicate substantial fluid motion. By com- 
parison, Figure 9 for attitude hold eight hours after the 
cessation of 1 rph rotation shows little evidence of fluid 
motion. 


From the foregoing we conclude that a simplified 
conductive model closely models attitude hold conditions 
and can serve as a useful real time indicator of the level 
of convective motion in the tanks. 
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CONVECTIVE BLOWDOWN 


The energy input required to maintain tank 
operating pressures of 900 psi increases sharply as low 
tank quantities are approached, and the dilute gas offers 
a diminished energy sink for the heater. Thus high heater 
temperatures may be encountered. An alternative operational 
mode is to take advantage of the pressurization energy and 
operate the tank as a compressed gas source. A constraint 
is that tank pressure be maintained at greater than 150 psi 
for proper flow to the fuel cells. One way of satisfying 
this requirement is to add sufficient energy to follow an 
isothermal path once past the apex of the two phase region, 
for example the -160° F isotherm in Figure 13. The energy 
required for total evacuation from 20% (71.3 lb) via such 
an isothermal blowdown is given in Figure 14 as a function 
of tank outflow rate. A portion of the needed heat is pro- 
vided by heat leak from the surrounding environment. In 
fact, the heat leak supplies all necessary energy for 
evacuation at an outflow rate of 0.5 lb/hr, and the heater 
duty cycles for 1.0 and 1.5 lb/hr flows are minimal ( < 1 0 % ) . 


CONCLUSIONS 


Two simplified models have been found to yield 
useful insights into the oxygen tank thermal behavior. The 
two fluid stratification model indicates that stratification 
of less than 1% of the fluid will lead to the observed 
pressure response times. The involvement of such a small 
fluid volume explains the insensitivity of the pressuriza- 
tion times to overall tank fluid motion and the rapid pres- 
sure recovery following fan operation. The conductive/ 
radiative heat transfer model furnishes a worst case heater 
temperature prediction and can serve as an indication of 
the extent of convective mixing. It is seen that the g- 
force exerted during attitude hold is insignificant in 
terms of enhanced heat transfer, but that ^3 rph passive 
thermal control has a noticeable effect. 
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FIGURE 5 - SCHEMATIC OF TWO FLUID MODEL 


195 



TANK PRESSURE (PSD 



FIGURE 6 • PRESSURE CYCLES BEFORE/AFTER FAN CYCLE 
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TX « TEMPERATURE OF NODE "X" 

CX = MASS X SPECIFIC HEAT OF NODE "X" IN BT0/°F s (VOLUMEI tpCp) 

GX = EFFECTIVE CONDUCTANCE BETWEEN ADJACENT NODES IN BTU/HR °F =kA/x 
G1 * EFFECTIVE RADIANCE BETWEEN NODES (1) AND (IS) IN BTU/HR °F 4 = Aea 


FIGURE 8 NETWORK DIAGRAM OF OXYGEN TANK THERMAL MODEL 
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TANK # 3 HEATER TEMP. (°F) 



FIGURE 9 - HEATER TEMPERATURE PROFILE DURING ATTITUDE HOLD 
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TANK = 3 HEATER TEMP, f F) 



FIGURE 10 HEATER TEMPERATURE PROFILE DURING ' 1 RPH SPACECRAFT ROTATION 
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FIGURE 11 - HEATER TEMPERATURE PROFILE DURING ' 3 RPH SPACECRAFT ROTATION 
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APOLLO 14 G E T. 

FIGURE 12 * HEATER TEMPERATURE PROFILE FOLLOWING SPACECRAFT ROTATIONAL MANEUVER 
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FIGURE 13 • OXYGEN SUPPLY PRESSURE VERSUS TANK QUANTITY 


203 


HEATER ENERGY INPUT IBTU/HR) 
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FIGURE 14 HEATER ENERGY INPUT FOR ISOTHERMAL BLOWDOWN FROM 20% QUANTITY 
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